Non-Fermi-liquid behavior due to short range order 
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An exactly soluble one-dimensional model of electrons interacting with order parameter flucta- 
tions associated with short-range order is considered. The energy and momentum dependence of 
the electronic self energy and spectral function are calculated and found to exhibit non- Fermi-liquid 
features similar to that seen for the two-dimensional Hubbard model: a pseudogap, shadow bands, 
anomalies in the self energy, and breakdown of the quasiparticle picture. Deviations from Fermi 
liquid behavior are largest close to the Fermi surface and as the correlation length increases. 

To appear in Physical Review B, Rapid Communications, November 1, 1996. 
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The question as to whe.ther the metallic properties 
of high- T„ superconductorsd and quasi-one-dimensional 
materialsH can be described by the Fermi liquid picture 
that works so well in conventional metals has recently re- 
ceived considerable attention. Many theoretical studies 
have been made of strongly correlated electron models 
such as the Hubbard and t — J m0 4fife§l-| The availability 
of high quality photoemission dataucliJ'Ll has recently fo- 
cusscd attention on the electron spectral weight function 
A(k, E) which is related to the probability of observing 
an electron with momentum k and energy E. 

A consistent picture is gradually emerging from studies 
of the two-dimensional HaJplaacd model using quantum 
Monte Carlo simulationstffl-ElO and studies jusing—the 
fluctuation exchange (FLEX) approximationEaOEjO. 
Some common features are observed as the temperature 
is lowered, the hole doping is decreased, or the Coulomb 
repulsion U is increased. It hasJapen suggested that these 
changes have a common originlid: they correspond to an 
increase in £, the correlation leagth associated with short 
range antiferromagnetic ordered. 

The common non-Fermi-liquid features observed are: 
(1) As £ increases a pseudogap develops, i.e., there is 
a suppi^siQBj-aL-the density of states near the Fermi 
energyfflliy'ElMEJ. (2) As £ increases peaks in the elec- 
tron spectral function become smaller and broadenU. 
Near the Fermi surface a breakdown of the quasiparticle 
picture may occur. (3) On the Fermi surface the real part 
of the self energVp£(fc, E) can develop a positive slope at 
the Fermi energyEd. (4) The magnitude of the imaginary 
part ofj-Lhc self energy has a local maximum at the Fermi 
energy!!-]. (5) "Shadow hands" exist due to the incipient 
antiferromagnetic orderliU, i.e., in addition to the peak in 
the spectral function at the energy E and momentum k 
there is a much smaller peak at k + Q where Q = (jr, w) 
is the wave vector associated with antiferromagnetic or- 
der (and the nesting vector of the Fermi surface for 
half filling). This means the spectral function A(k,E) 
has peaks for \k\ > kp but E < Ep. It is estimated 



that such shadow bands are observable when the corpa- 
lation length is larger than a couple of lattice spacingsll3. 
The corresponding photoemission peaks, have been ob- 
served in metallic Bi2Sr2CaCu20s+M'" and insulating 
Sr2CuC>2Cl2u. (6) The seif energy associated with the 
shadow bands is singulania. As £ increases the imagi- 
nary part of the self energy is much larger on tho-shadow 
Fermi surface than on the regular Fermi surfaceO. 

The purpose of this Rapid Communication is to point 
out that there is a simple one-dimensional model which 
also has many of the features listed above. It is easy to 
study because, it has an exact analytic solution, found 
by Sadovskiilia This solution can be used to test differ- 
ent approximation schemes and methods of analytic con- 
tinuation used on two-dimensional models. Hopefully 
the results presented here will also provide physical in- 
sight into how short range order produces non-Fermi- 
liquid behavior in two dimensions. The model is directly 
relevant ta_the related issues for quasi-one-dimensional 
materialscjO. 

The model consists of left and right moving elec- 
trons, with Fermi velocity Vp, interacting with a static 
backscattering random potential A(z), with Hamiltonian 



H = / dz& 



d 

WF<J3 ~dz + ^ z ^ a+ + A(z)* cr - 



* (1) 



where a% and a± = ± are Pauli matrices. The 
upper and lower components of the spinor ^ (z) are left- 
moving, up-spin and right-moving, down-spin electrons, 
respectively. The other electrons are described by a sim- 
ilar Hamiltonian. (The form of the Hamiltonian is moti- 
vated by considering backward scattering in the Hubbard 
model). The random potential has zero mean and finite 
range Gaussian correlations given by: 



<A(z)A(zT}=^ 2 exp(-|z-z'|/0- 



(2) 



where ij> is the rms fluctuation in the potential at a given 
point and £ is the correlation length, ip defines an energy 
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scale and a length scale £o = vp/ip. It will be seen below 
that the ratio of the correlation length to £o determines to 
what extent the short-range order causes deviation from 
Fermi liquid behavior. For the commensurate case of a 
half-filled band A(z) is real. 

An alternative interpretation of the Hamiltonian is 
that it represents electrons interacting with spin fluctu- 
ations, S(x) = A(x)/tp, with susceptibility 



(q±2fc F )2+e- 



(3) 



and ip then defines the strength of the coupling of the 
electrons to the spin fluctuations. In this interpreta- 
tion we are assuming that the energy scale associated 
with the spin fluctuations is so much smaller than the 
electronic energy scale that the spin fluctations can be 
treated as static. Deisz, Hess, and SereneEj noted that 
for the two-dimensional Hubbard model at half filling 
anomalies in the self energy were always accompanied by 
a spin fluctuation T-matrix T aa {q,bj n ) that is strongly 
peaked near the wave vector q = Q and Matsubura fre- 
quency uj n = 0. In this regime there is a natural map- 
ping onto the model considered here. At half-filling the 
Fermi surface is square and has perfect nesting so is some- 
what "one-dimensional" . The one-dimensional momen- 
tum then corresponds to that along the (n, tt) direction 
in the two dimensional model, ip 2 corresponds to the 
weight of the peak iaJT acT (q, 0) within of Q (compare 
equation (5) in Ref.E-3). 

Sadovskii considered the model defined by (|l|) and (||) 
and found an exact solution by summing, all the dia- 
grams generated by perturbation theory!!! The elec- 
tronic Greens function is a continued fraction 



G(k,E) = 



1 



E - do - 



(4) 



E-ai - 



where for right moving electrons 



1p C 2 



E-a 2 



ai = (-1) v F k + sign(fc) — 



(5) 



and Ci = I for the commensurate case of a half-filled 
band and for the incommensurate case q — 1/2 for I 
even and ci — [l + l)/2 for I odd. As far as we are 
aware this is the only non-trivial electronic model for 
which there is an exact analytic solution for the con- 
tinued fraction representation of a correlation function. 
Since the Lanczos method of exact diagonalization pro- 
duces such a continued fraction representation of spectral 
functions^ this model could provide insights into trends 
in the continued fraction coefficients and possible termi- 
nation proceduresE3. In the limit £ — ► oo a perturbation 
expansion for G(k, E) is divergent but Borel summableEI 
We evaluated the continued fraction (|) numerically us- 
ing the modified Lentz's methodca. All the results shown 



here are for the incommensurate case. Qualitatively sim- 
ilar results are obtained for the commensurate case. 

Generally we find that as the correlation length in- 
creases and the Fermi surface is approached the prop- 
erties of the self energy and spectral function deviate 
significantly from the quasiparticle picture of Fermi liq- 
uid theory. Many of the anomalous features we see are 
similar to those found for the two-dimensional Hubbard 
model at half-JiJling treated in the fluctuation-exchange 
approximatiorO. Figure 1 shows how the spectral func- 
tion A(k, E) broadens significantly as the electron mo- 
mentum approaches the Fermi momentum. This happens 
even for correlation lengths of the order of £o = vp/ip. 
This means that for strong coupling, i.e., ip of the order 
of the band width, the correlation length £ can be of the 
order of a lattice constant. Shadow bands are present 
and become larger as £ increases and \k\ decreases. 

Figure 2 shows how the quasiparticle picture breaks 
down as the correlation length £ increases. As £ increases 
from 0.2£o to 100£o the spectral function on the Fermi 
surface evolves from a single narrow peak to two broad 
peaks which are the precursors of conduction and va- 
lence bands associated with long range spin-density-wave 
order. For similar reasons as the correlation length in- 
creases a pseudogap develops in the total density of states 
(Figure 2 inset). It should be pointed out tiiat the two 
peak structure is not seen in FLEX resultsfi-J or quan- 
tum Monte carlo for weak coupling on large latticesB. 
Whether this absence of the two peaks is a real prop- 
erty of the 2D Hubbard model or a result of the FLEX 
approximation or not being able to go to low enough 
temperatures is not clear. 

The self-energy has anomalous features similax.to those 
found for the 2D Hubbard model at half-fillingEl At the 
Fermi energy the real part of the self energy has a positive 
slope (Figure 3). The magnitude of the imaginary part of 
the self energy, which in a quasiparticle picture is related 
to the quasiparticle lifetime and develops a maximum 
at the Fermi energy (Figure 3). In contrast, in a Fermi 
liquid ReE has a negative slope and |ImI]| has a minimum 
(which goes to zero as the temperature goes to zero) at 
the Fermi energy. The opposite behavior observed here 
means that although for our model there is a pole in the 
spectral function at E = it is not possible to define a 
quasiparticle solution there. 

The self-energy on the Fermi surface S(0, E) develops 
a singularity at the Fermi energy E = as the correlation 
length increases. Figure 4 shows that T,(k,E) develops 
similar singular behavior near E ~ —kvp which corre- 
sponds to the position of the shadow band. In general, 
the magnitude of the imaginary part of the self energy 
(which is related to the scattering rate) is much larger 
and more singular near E ~ —kvF than E ~ +kvF- 
The scattering rate for the shadow band increases with 
£, as is observed for the 2D Hubbard model as the dop- 
ing or temperature is decreasedEj. Only for very large 
correlation lengths does the shadow band feature corre- 
spond to a pole of the spectral function, i.e., an energy 
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Ek which satisfies E^ — kvp — Re£(fc,£^) — 0. Con- 
sequently, the shadow band peak is not a replica of the 
regular quasiparticle peak—These results can be com- 
pared to those of Chubukovc3 for a two-dimensional spin- 
fluctuation model. 

Figure 5 shows the spectral function calculated for 
this model using second order perturbation theory. This 
corresponds to termination of the continued fraction at 
1 = 1. Such a perturbative form for the Greens func- 
tion has been used in calculations concerning the role 
of order parameter fluctuations in quasi-one-dimensional 
materials£-3. However, the discrepancy with the exact re- 
sults shown in Figure 1 shows this is quite unreliable for 
\E\ < ip when \k — kp\ < tp/vp and £ > Vp/ip. In particu- 
lar perturbation theory greatly underestimates the width 
of the spectral function. The inset of Figure 5 shows how 
perturbation theory gives unreliable results for the total 
density of states if £ > £o ■ 

In conclusion, the simple model considered here 
has many features similar to those seen for the two- 
dimensional Hubbard model and provides insight into 
how short range order can produce non-Fermi-liquid be- 
havior. Of particular interest is the following, (a) The 
ratio of the correlation length to £o — vp/ip determines 
deviations from Fermi liquid behavior, (b) Non-Fermi- 
liquid behavior occurs even when the correlation length 
is sufficiently short that there is only a weak pseudogap. 
(c) The "shadow band" is associated with singularities 
in the self energy and only corresponds to poles in the 
spectral function for very large correlation lengths, (d) 
Perturbation theory gives unreliable results in the non- 
Fermi-liquid regime. 

This work was stimulated by discussions with J. R. 
Schrieffer. We thank H. Castella, J. Deisz, D. Hess, M. 
Steiner, and J. Voit for very helpful discussions. This 
work was supported by the Australian Research Council. 

Note added in proof. Recent angle-resolved photoe- 
mission measurements [H. Ding et ai, Nature 382, 51 
(1996); A. G. Loeser et ai, Science 273, 325 (1996)] 
have measured a pseudogap in the normal state of 
Bi^SriCaCu-iO^+x. The observed opening up of the 
pseudogap with decreasing doping and temperature are 
qualitatively similar to the how the pseudogap considered 
here opens up with increasing correlation length (Fig. 2 
inset). 
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FIG. 1. Shadow bands and the breakdown of the quasi- 
particle picture near the Fermi surface. The energy depen- 
dence of the electron spectral function A(k,E) is shown for 
three different momenta k. The correlation length associ- 
ated with the short range order is equal to 5£o = 5vf/^> 
where vf is the Fermi velocity and tp is a measure of the 
strength of the coupling of the electrons to the order param- 
eter fluctuations. Note the presence of "shadow bands", i.e., 
small peaks in A(k, E) with E < Ef = and k > kF = 0. 
The spectral function becomes significantly broader as the 
Fermi surface is approached. Opposite behavior occurs in a 
Fermi liquid: the quasiparticles are better defined close to the 
Fermi surface. These results can be compared to those for the 
two- dimensional Hubbard model at half filling (see Figure 1 
in RefH and Figure 2 in Ref.H). 
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FIG. 2. Breakdown of the quasiparticle picture with in- 
creasing correlation length. The energy dependence of the 
spectral function at the Fermi momentum is shown for sev- 
eral correlation lengths. The spectral function broadens con- 
siderably and evolves into two bands as the correlation length 
increases. These results can be compared to those for the 2D 
Hubbard model at half-filling as the temperature is lowered 
(see Figure 1 in Ref.E-3). The inset shows how the pseudo- 
gap in the total density of states opens up as the correlation 
length increases. Note that although for f ~ £o the pseudo- 
gap is rather weak deviations from Fermi liquid behavior still 
occur. 
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FIG. 3. Anomalous behaviour of the self energy. The real 
part of the self energy at the Fermi momentum, ReE(0,iS), 
has a positive slope at the Fermi energy (E = 0). This means 
that a quasiparticle weight cannot be defined. As the cor- 
relation length increases a maximum at E — develops in 
the magnitude of the imaginary part of the self energy atrfelpe 
Fermi momentum |lm£(0, E)\. (Compare Figure 2 in Rem). 



4 



-1 





E/V 



1 



FIG. 4. Singular behavior of the shadow band self energy. 
Away from the Fermi surface the self energy near E ~ kvF has 
Fermi liquid properties. In contrast, near E ~ —kvF which 
is associated with the shadow band the self energy becomes 
increasingly singular as the correlation length increases. In- 
tersections of the ReE curves with the dot-dashed straight 
line correspond to poles of the spectral function. Hence, the 
shadow band is only associated with a pole for very large 
correlation lengths. 
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FIG. 5. Second order perturbation theory is unreliable near 
the Fermi surface. The spectral function is calculated for the 
same parameter values as in Figure 1. (Note the vertical 
scale is three times larger here). The discrepancy with the 
exact results is large for \k\ < l/£o and becomes larger as £ 
increases. The inset shows the total density of states for the 
same parameter values as in Figure 2. (Note the vertical scale 
is two times larger here). 
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